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Abstract. We present BVRI photometry and spec- 
trophotometry of GRB 990510 obtained with the ESO 
VLT/ Antu telescope during the late decline phase. Be- 
tween days 8 and 29 after the burst, the afterglow faded 
from R = 24.2 to ~ 26.4. The spectral flux distribution 
and the light curve support the interpretation of the after- 
glow as synchrotron emission from a jet. The light curve is 
consistent with the optical transient alone but an under- 
lying SN with maximum brightness R > 27 .4 or a galaxy 
with R > 27.6 (3-er upper limits) cannot be ruled out. To 
a 5-<7 detection threshold of R = 26.1, no galaxy is found 
within 6" of the transient. A very blue V ~ 24.5 extended 
object which may qualify as a starburst galaxy is located 
12" SE, but at unknown redshift. 

Key words: Gamma-rays: bursts - optical radiation - 
Acceleration of particles - Galaxies: general - Galaxies: 
starburst - Cosmology: miscellaneous 



1. Introduction 

Gamma-ray bursts (GRBs) have their origin in highly dy- 
namical processes which result in rclativistic blast waves 
(e.g. Piran 1999). The currently best candidates are merg- 
ers of compact stars or the collapse of massive stars. The 
blast wave and its interaction with the surrounding in- 
terstellar or circumstellar matter creates an optical after- 
glow which has been observed for eleven GRBs thus far. 
For seven of these, the host galaxies have been found and 
for five the redshift could be measured (Hogg & Fruchter 
1999). While the final cataclysms of massive stars occur 
within or near the host galaxies, less massive systems like 
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binary neutron stars can be expelled from their hosts at 
large velocities and suffer the merging event at distances 
from the hosts exceeding a Mpc (Fryer et al. 1999). 

GRB 990510 was detected by BATSE, Ulysses, and 
BcppoSAX as a very bright burst lasting 68 s (see Wi- 
jers et al. 1999 for the burst profile and a first summary). 
The optical afterglow was discovered by Vreeswijk et al. 
(1999a) and subsequently followed by numerous observers. 
An early spectrum located the optical transient (OT) at a 
redshift z > 1.62 (Vreeswijk et al. 1999b). GRB 990510 is 
the first burst to show a clearly defined achromatic break 
in the optical light curve (Harrison et al. 1999, Israel et 
al. 1999, Stanek et al. 1999) which was readily interpreted 
as firm evidence for beaming and a total energy release 
substantially less than the isotropic value of 3 x 10 53 ergs. 
In this Letter, we report observations with the ESO VLT 
of the late decline of the optical afterglow at two epochs. 

2. Observations 

GRB 990510 was observed with ESO's Very Large Tele- 
scope (VLT) Unit 1 Antu equipped with the Focal Reducer 
Low Dispersion Spectrograph (FORS1) between May 14 
and 18, 1999, and again between June 8 and 11, 1999 
(3.8-8.0 and 28.7-31.6 heliocentric days after burst, here- 
after HDAB). We performed photometry in the Bessel 
B,V,R, and the Gunn I bands, supplemented by long-slit 
spectrophotometry. The scale was 0.2" /pixel in all ob- 
servations except for the B, R images on June 8, 1999 
which were taken in the high-resolution mode of FORS1 
with 0.1" /pixel. The data were reduced with standard 
ESO-MIDAS procedures. Photometry was derived rela- 
tive to the comparison stars given by Pietrzynski & Udal- 
ski (1999), with stars A and B near the optical transient 
(Fig. 3) serving as secondary standards (Table 1). Due 
to wind buffeting, the point spread function (PSF) was 
slightly non-circular in some images, so magnitudes were 
determined by a special PSF-template routine. Table 1 
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provides a summary of the observations with times given 
in heliocentric Julian days (HJD) and HDAB. 1-a internal 
statistical errors arc quoted for the photometric magni- 
tudes. 

Long-slit spectrophotometry was performed on May 
14 and 16, 1999. The slit was oriented at a position angle 
— 15° to include objects C south and A north of the OT 
(Fig. 3) and, coincidcntally, also two M stars 12" north 
of the OT (R = 22.1,23.5, outside of Fig. 3). The late 
spectroscopy of the OT is reported elsewhere (Vreeswijk 
et al. 1999). 

The seeing encountered during the observations was 
never better than 0.8" (Table 2). This resulted in 5- 
a detection thresholds of B = 26.6 on June 8, V = 
25.6 on May 18, R = 26.1 on June 8, and / = 25.0 
on May 14 and June 11. 



3. Results 

3.1. Spectrophotometry 

Our slit spectra of objects A and C reveal their stellar 
nature. Object A is a dK6 ± 1 star of the disk population. 
Its brightness (Table 1) places it at a distance of d ~ 2 kpc, 
outside the galactic dust layer. The spectral type implies 
an intrinsic V - 1= 1.40 ± 0.13 and £ B -v = 0.14 ± 0.09, 
consistent with Eb-v= 0.20 ± 0.03 quoted by Stanek et 
al. (1999). Object C, initially suspected to be the host 
galaxy, is pointlikc. Its spectrum is included in Fig. 1. 
The probable presence of Mgl A5167 suggests that it is an 
early dK star at a distance of ~ 15 kpc. The photometry of 
objects B and D suggests they are stars of spectral types 
dK5 and dMO, respectively (Table 1). 

Figure 1 also shows the extinction-corrected spec- 
trophotometric flux distribution of the OT on HDAB 3.9 
when the OT had R = 21.98, corrected with £ B -v = 0.20 
(histogram; binned in intervals of 200 A). The solid dots 
denote the photometric result. For the interval of 4900- 
9000A, the flux distribution can be described by a power 
law F v oc v'P oc X fi with f3 = 0.55 ± 0.10. For higher fre- 
quencies, the spectral flux distribution steepens, as noted 
already by Stanek et al. (1999). The steepening may be 

Table 1. Bessel photometry of stars A, B, C, and D in 
Fig. 3. The errors quoted refer to the internal statistical 
uncertainties only. The primary standard is star No. 3 
of Pictrzynski & Udalski (1999) with B = 17.88, V = 
17.01, R = 16.50, and / = 16.09. 
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Fig. 1. Spectral flux distribution of the OT of 
GRB 990510, de-reddened with E B -v = 0.20. The solid 
points refer to the mean magnitudes of the OT on HDAB 
3.9, the histogram presents our spectrophotometry aver- 
aged over 200A bins. The dashed line is the power-law 
fit to A > 4900A with slope f3 = 0.55 for a power law 
F v oc oc X 13 . Shown below is the spectrum of the 

V = 22.6 early K-star C, 2" south of the OT, shifted 
downward by a factor of two. 

intrinsic to the source or indicate additional absorption 
outside our Galaxy. 

3.2. Light curve 

Our BVRI photometry provides colour information late 
in the afterglow which is complementary to the ample 
information gathered earlier by other observers. Around 
HDAB 1.5, the light curve was found to change from an 
early power law F = k t~ a with a = 0.8 to a steeper power 
law with index 2.2 (Harrison et al. 1999, Israel et al. 1999, 
Stanek et al. 1999, and references therein). 

In the HDAB 0.6-1.1 interval, i.e. before the break, the 
OT had colours B-V = 0.57 ±0.02, V-R= 0.41 ±0.01, 
and R—I = 0.47±0.01. The early spectral flux distribution 
was possibly slightly bluer, as judged from the Mt. Stromlo 
data (Harrison et al. 1999) which yield V-R = 0.31 ±0.03 
at HDAB 0.15 and V - R = 0.36 ± 0.05 for HDAB 0.4. 
Throughout the further evolution, however, there is no 
evidence for a significant change in colours: (i) in the 
HDAB 1.6-2.0 interval are B — V = 0.61 ± 0.07, V — R = 
0.42 ± 0.03, and R - I = 0.40 ± 0.02, consistent with the 
colours before the break; (ii) on HDAB 3.85, well after 
'the brak, our BRI photometry yields B-R= 0.98 ± 0.07 
■and R — I = 0.49 ± 0.06; and (iii) the latest colour in- 
formation available, for HDAB 5.7-8.0, yields a mean 

V — R = 0.37 ± 0.03 (all colours derived from magnitudes 
corrected for long-term trends) . These results demonstrate 
that the decay of the afterglow is achromatic in VRI at 
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least for the interval HDAB 0.6-8.0 and in BVRI for 
HDAB 0.8-3.9. We feel justified, therefore, in converting 
all measurements to equivalent i?-magnitudes and fitting 
a single grand total light curve (Fig. 2. upper panel). The 
final mean colours used in this conversion were determined 
iteratively and are given below. 

The errors attached to the individual data points are 
the statistical errors derived from the noise in the CCD 
images. Systematic zero point errors are typically quoted 
as 0.02-0.03 mag, consistent with the internal scatter of 
~ 0.02 mag and the lack of obvious time variability among 
the 101 (of a total of 163) data points in the HDAB 0.6-1.1 
interval. While short term variability is certainly small, a 
systematic search for such variability among all photomet- 
ric data is still pending. In order to see if systematic zero 
point errors affect the fit, we also analysed a data set for 
which 0.03 mag was quadratically added to the statistical 
errors before fitting the light curve. 

Burst models (Piran 1999, Rhoads 1999, Sari et al. 
1999) explain the afterglow as synchrotron emission of 
shock-accelerated electrons injected into an expanding me- 
dium with a power law spectrum E~~ p . In jet models, the 
time dependent spectral flux at frequencies below the cool- 
ing break varies as F v (t) oc v^ 13 t~ a , with (3 = (p — l)/2 
independent of time and a = 3(p - l)/4 or a = p, de- 
pending on, respectively, whether the opening angle of 
the relativistically beamed radiation d ~ 1/7 < 6 early 
in the expansion, or 1? > 6 at later time when the jet has 
been slowed down (8 = opening angle of the jet, 7 = bulk 
Lorentz factor). We choose a function 

F(t) = (Ff" + Fa - ")- 1 /" with Ft = hr a % n > 0(1) 

to describe the transition between the early and late power 
laws Fi and F2, where F\ — F2 at the transition time 
t = t*. Eq. (1) was also employed by Rhoads (1999) to 
parameterise his numerical models and is a more general 
form of the expressions used by Israel et al. (1999), Stanek 
et al. (1999), and by Harrison et al. (1999), who assumed 
n = 1 and n ~ 1.5, respectively. Of the five free param- 
eters in Eq. (1) (fci,2, cci,2, n), the exponent n provides a 
measure of the relative width and the smoothness of the 
transition from F\ to F2. The extrapolation of the OT 
brightness from HDAB< 8.0 to late times will also de- 
pend on the choice of n. 

We first fit the light curve of the OT using all of the 
available data up to HDAB 8.0. For the data sets with 
zero (0.03 mag) systematic errors added, we find ct\ = 
0.79±0.04 (0.81±0.05), a 2 = 2.41±0.16 (2.31±0.16), n = 
0.87±0.22 (1.09±0.34), and t* = 1.75±0.23 (1.57±0.19) 
days. The quoted uncertainties are 1-ct errors. The fit- 
ted x 2 /^ = 333/165 d.o.f = 2.02 shows that the fit to 
the unaltered data is not perfect, but is acceptable if sys- 
tematic errors of 0.03 mag are considered (x 2 /^ = 0.83). 
The observed n ~ 0.9 — 1.1 differs significantly from the 
value n = 0.4 found by Rhoads (1999) when fitting his 



numerical results, which predict a rather smooth tran- 
sition. The more abrupt break may be due to the su- 
perposition of several effects, as the slowing of the jet 
and its sideways spread (Harrison et al. 1999). By fit- 
ting the shape of the mean light curve to the fluxes in 
each photometric band we obtain the magnitude normal- 
isations corresponding to the quantities fc, in Eq. 1: B\ = 
20.550 ±0.012, Vi = 19.961 ±0.003, R x = 19.565 ± 0.003, 
and I\ = 19.075 ± 0.003. The corresponding colours are 
consistent with those for HDAB 0.6-1.1. The lack of sub- 
stantial colour variations can be seen from the small scat- 
ter in the residuals of the best fit (Fig. 2). 

For comparison, for n fixed at 1.0 we obtain ot\ = 
0.80 ± 0.02, a 2 = 2.33 ± 0.04, t* = 1.63 ± 0.08 days and 
X 2 /^ = 2.0 for a fit to the unaltered data points. Since 
the errors of the a's and of n are correlated, fixing n re- 
duces the errors in the other parameters artificially. Our 
fit to the light curve with n ^ 1 allows us to predict the 
brightness of the OT at later times with some confidence 
that uncertainties in the shape of the light curve have been 
accounted for as far as possible. We predict an OT mag- 
nitude of R = 26.92-^j? on June 8 (HDAB 28.8), when 
our BR photometry was performed. 

Fruchter et al. (1999) have reported detections with the 
HST-STIS open CCD on June 8.1 and June 17.9 (HDAB 
28.7 and 38.5). Assuming an F„oc spectrum, they 

obtained V = 27.0±0.2 and V = 27.8±0.3, corresponding 
to R = 26.6 ± 0.2 and R = 27 A ± 0.3, respectively. Given 
the red sensitivity of the STIS-CCD, these R magnitudes 
are quite robust. We conclude that our June 8 result and 
the first of the HST detections are consistent with each 
other. 

Assuming that all photometry including the late data 
points is that of the decaying OT, we obtain nearly identi- 
cal light curve parameters: R\ = 19.539 ± 0.063 (19.555 ± 
0.082), ai = 0.82 ± 0.03(0.84 ± 0.03); a 2 = 2.23 ± 
0.07 (2.18 ±0.07); n = 1.17 ± 0.17 (1.43 ± 0.26); U = 
1.51 ± 0.09 (1.41 ± 0.08); and X 2 /v = 1.87 (0.81). Thus, 
all of the available photometry is consistent with the light 
curve of a pure OT. 

The observed spectral index (3 ~ 0.55 and the values 
of «i and a 2 are in excellent agreement with the theoret- 
ical predictions if the cooling break stays at wavelengths 
shorter than the visual band until HDAB 8.0. The slope 
of the electron injection spectrum then is p = 2.1 ± 0.1 in 
agreement with the conclusions drawn by Harrison et al. 
(1999) and Stanek et al. (1999). 

Figures 3a-c depict the late decay of the OT from an 
equivalent i?-magnitudc of 21.95 (HDAB 3.8) over 23.75 
(HDAB 8.0) to 26.4 (HDAB 28.8) (Table 2). Down to 
the 5-er thresholds for a detection at an arbitrary posi- 
tion, B ~ 26.6, V ~ 25.6, R ~ 26.1 and I ~ 25.0 our 
images show no other object within 6" of the OT, be- 
sides the stars A, B, C, and D. The faint emission on 
June 8 at the OT position is more easily seen when the 
high- resolution image (0.1" pixels) is slightly smoothed. 



K. Beuermann et al.: VLT observations of GRB 990510 



QJ 

■+-> 

. i — i 

C 

w 

CO 



18 



20 - 



22 



24 



26 - 



28 



-0.5 
0.0 



10" 



GRB990510 



30 d 40 d 50 d 60 d 








R SN" 


i 

= 27.4 ~ 






/ 








\ 


R gal -27.6 : 






-* -N 




































LATE 


OT - 


\ V. 











Residuals for Early OT Fit 



_J I I I I I I 1 I I L I I LJ_ 



10 u 10 1 
Days After Burst 



10" 



Fig. 2. Top: Light curve of the OT of GBR990510, with all published measurements converted to the i?-band. The 
solid curve indicates the fit to HDAB< 8.0 data. Insert: Fits to all available data assuming the OT alone (Late OT, 
dotted), the OT with a 3-cr upper limit flux to a supernova (long dashes) or a host galaxy (short dashes), x HST, o 
late VLT data. Bottom: residuals relative to a pure OT light curve fit to the HDAB < 8 data points. The 3-cr bounds 
(dashed lines) show the uncertainties in the extrapolation to late times. The dotted line indicates the pure OT fit 
including the three late data points. 



Table 2. Log of our observations of GRB 990510 with the ESO VLT Antu (UT1) and the Focal Reducer Spectrograph 
FORS1. HDAB is the time of the centre of the exposure in heliocentric days after the burst which occurred on HJD 
2451308.870. \-a statistical errors are quoted. 



Date of 


HJD 


HDAB 


Observ. 


Seeing 


Filter or 


Expos. 


Flux rel. to star A 


Magnitude 


Equiv. R-mag 


1999 


-2451300 


(days) 




(arcsec) 


Grating 


(sec) 


(%) 






May 14 


12.707 


3.837 


Phot. 


0.8 


R 


300 


7.67 ±0.11 


21.95 ±0.02 




May 14 


12.716 


3.846 


Phot. 


0.9 


I 


1200 


6.24 ±0.07 


21.45 ±0.01 


21.93 ±0.01 


May 14 


12.730 


3.860 


Phot. 


1.1 


B 


100 


23.01 ±0.96 


22.93 ±0.05 


21.95 ±0.05 


May 14 


12.763 


3.893 


Spectr. 


0.8 


G300V 


1800 










12.789 


3.919 


Spectr. 


0.8 


G300I/OG590 


1800 








May 16 


14.731 


5.861 


Phot. 


1.2 


V 


600 


4.41 ±0.14 


23.39 ± 0.04 


22.99 ± 0.04 


May 16 


14.746 


5.876 


Spectr. 


1.5 


G150I 


1800 








May 18 


16.867 


7.997 


Phot. 


1.0 


V 


1200 


2.24 ±0.03 


24.16 ±0.06 


23.75 ±0.06 


June 8 


37.628 


28.758 


Phot. 


0.8 


B 


1200 


0.57 ±0.20 


> 26.7 


> 25.7 


June 8 


37.656 


28.786 


Phot. 


0.9 
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2400 


0.14 ± 0.04 


26.34 ±0.27 




June 11 


40.485 


31.615 


Phot. 


0.9 


I 


3000 


0.13 ±0.09 


> 25.1 


> 25.6 
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Fig. 3. (a) - (c), from left: Deep I, V, and R images of the field of GRB 990510 taken on May 14, May 18, and June 
8 when the OT was at I = 21.46, V — 24.16, and R ~ 26.7. (d) Smoothed i?-image of June 8 with star C subtracted 
to better show possible nearby faint sources. The circles at the positions of A, C and the OT have radii of 1 FWHM 
of the PSF. Object E is a 24 mag possible starburst galaxy, 12" SE of the OT. 



Figure 3d is an enhanced version of 3c, smoothed with 
a Gaussian of 0.3" standard deviation and with star C 
PSF-subtracted. While the emission appears to be slightly 
offset from the OT position by 0.4" and may signify an 
underlying object, the shift is within the statistical uncer- 
tainty expected for such a faint source. 



3.3. The quest for the host galaxy of GRB 990510 

Although the light curve analysis does not require another 
light source in addition to the OT, such contribution can 
not be excluded. One possibility is that an underlying 
SN-event contributes to the June 1999 light level. Tak- 
ing the i?-band light curve of SN1998bw (Galama et al. 
1998, Iwamoto et al. 1998) adjusted to z = 1.62 (Bloom 
et al. 1999) one expects a peak magnitude of R = 26.8 
(reddened 27.3), just consistent with the 3-cr upper limit 
of R — 27 .4 of the fit [fitted flux relative to star A with R 
= 19.17 is (1.0 ±1.4) x 10~ 4 ]. Forcing the fit to use a peak 
value of 27.4 results in an increase in ct2 from 2.37 to 2.71 
and a correlated decrease in n from 0.94 to 0.60. Alterna- 
tively, the best-fit contribution by an underlying galaxy, 
which adds a constant to the late-time light curve, yields a 
3-cr upper limit of 27.6 (27.1 de-reddened) [flux relative to 
star A (2.0±1.7) x 10~ 4 ]. The forced fit with the 3-cr upper 
limit requires cx-i and n changes nearly identical to those of 
the corresponding SN fit. We note that many burst models 
do not require a SN event and that a host fainter than 27 
mag presents no problem with the galaxy population re- 
sponsible for GRBs (Hogg & Fruchter 1999). Fig. 2 shows 
the fits obtained for the OT alone (early data or all data) 
and the forced fits with the OT and the 3-cr upper-limit 
SN and host contributions added: there is presently no 
way to distinguish between the three models, but the fact 
that the early and late OT light curve fits are so similar 
strongly suggests that we have only seen the OT. 



For z — 1.62 (Galama et al. 1999) and a standard cos- 
mology with H Q = 70kms~ 1 Mpc~ 1 and £1 = 0.3, the 
distance modulus is 45.2. The de-reddened host magni- 
tude R > 27.1 then translates to a restframe Mr > —18.4 
for a starburst spectrum similar to NGC4449 (Bruzual 
& Chariot 1993) and to Mr > —17 for a pure starburst 
spectrum. Hence, the host is intrinsically faint if at z= 1.62. 
While one might be willing to place the burst at a larger z, 
there is a probable limit from the expected Lya-forest or 
Lyman-continuum absorption. The lack of a pronounced 
depression at wavelengths longer than 3600A (Fig. 1) lim- 
its the acceptable redshift to z Si 2.0 for Lya-forest and 
to z Si 2.9 for Lyman-continuum absorption. 

The only brighter object possibly related to the OT 
is the one denoted E in Fig. 3d, 12" to the SE. With a 
diameter of ~ 1" and its decidedly blue colour it may 
qualify as a starburst galaxy. Its de-reddened magnitudes 
and colours are well within the range found for the hosts of 
other GRBs and OTs, V ~ 23.8, B-V ~ 0.1, V-R ~ 0.3, 
R—I ~ 0.4, assuming Eb-v = 0.20. The observed spectral 
flux distribution is typical of the UV spectrum of a star- 
burst or irregular galaxy (e.g. Bruzual & Chariot 1993). 
If this is actually a galaxy at z = 1.62, it would be sepa- 
rated by at least 100 kpc from the OT. This would not be 
a problem if the GRB progenitor were a low-mass system 
ejected from this galaxy (Fryer et al. 1999). The absence of 
an associated SN or host galaxy at the OT position would 
support such a scenario, while the suspected connection 
of long gamma-ray bursts with massive progenitors does 
not. 

4. Conclusions 

The recent observations of the optical transient of 
GRB 990510 strongly support the synchrotron model for 
gamma-ray bursts. The pre-break light curve and the 
A > 4900 A optical flux distribution arc consistent with 
the adiabatic cooling of relativistic electrons with an E~ p 
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injection spectrum with p ~ 2.1. The steep late decline is 
a firm indicator of the presence of a slowing jet (Meszaros 
& Rees 1999, Rhoads 1999, Sari, et al. 1999). 

The light level on June 8 is quite consistent with being 
solely due to the transient, but contributions by either 
a SN with R > 27 A or a host with R > 27.6 (3-cr up- 
per limits) cannot be ruled out. If located at z = 1.62, 
the SN could have been nearly as bright as expected from 
appropriately scaling SN1998bw, while a starburst galaxy 
as the host would have Mr^ —17. Proving the presence 
or absence of a host galaxy at the OT position or iden- 
tifying fainter nearby candidates requires additional deep 
exposures. If ejection of the progenitor of GRB 990510 is 
considered a possibility, the blue extended object 12" east 
with de-reddened R — 23.6 could be a host candidate. 
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